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Abstract 

Injuries to articular cartilage result in the development of lesions that form on the surface of the cartilage. 
Such lesions are associated with articular cartilage degeneration and osteoarthritis. The typical injury 
response often causes collateral damage, primarily an effect of inflammation, which results in the spread of 
lesions beyond the region where the initial injury occurs. We present a minimal mathematical model based 
on known mechanisms to investigate the spread and abatement of such lesions. In particular we represent 
the "balancing act" between pro-inflammatory and anti-inflammatory cytokines that is hypothesized to be 
a principal mechanism in the expansion properties of cartilage damage during the typical injury response. 
We present preliminary results of in vitro studies that confirm the anti-inflammatory activities of the 
cytokine erythropoietin (EPO). We assume that the diffusion of cytokines determine the spatial behavior 
of injury response and lesion expansion so that a reaction diffusion system involving chemical species 
and chondrocyte cell state population densities is a natural way to represent cartilage injury response. 
We present computational results using the mathematical model showing that our representation is 
successful in capturing much of the interesting spatial behavior of injury associated lesion development and 
abatement in articular cartilage. Further, we discuss the use of this model to study the possibility of using 
EPO as a therapy for reducing the amount of inflammation induced collateral damage to cartilage during 
the typical injury response. The mathematical model presented herein suggests that not only are anti- 
inflammatory cytokines, such as EPO necessary to prevent chondrocytes signaled by pro-inflammatory 
cytokines from entering apoptosis, they may also influence how chondrocytes respond to signaling by 
pro-inflammatory cytokines. 
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Introduction 

Articular cartilage is composed of cells known as chondrocytes. Mechanical stress and injury is known to 
kill (via necrosis) chondrocytes and results in the formation of lesions on the cartilage surface [l}[4] . The 
limited capacity of chondrocytes to self-repair, together with certain aspects of the typical injury response 
such as inflammation, can cause the spread of lesions and development of osteoarthritis. Recent research 
[5] suggests that inflammatory cytokines such as tumor necrosis factor a (TNF-a) play a significant role 
in causing the spread of cartilage lesions. Anti-inflammatory cytokines such as erythropoietin (EPO) 
play an antagonistic role to TNF-a, limiting the area over which a lesion can spread by counteracting 
some of the effects of inflammation [sj. Moreover, the authors in [5] describe the potential use of EPO 
as a therapy for cartilage injury and lesion abatement 

The potential for EPO to suppress catabolic and inflammatory responses of chondrocytes to various 
alarmins, i.e. molecules that trigger the innate immune response, was explored in a cell culture model 
(Figure [1]). Bovine articular chondrocytes in serum-free culture were pre-treated with EPO and then 
exposed to alarmins in the form of cell lysates or the 29 kDA fibronectin fragment (Fn-f). Intcrleukin-1/3 
(IL-1/3) was used as a control. Cell lysates were made by repeated freeze-thawing of excess cells from 
the same culture. Culture media were collected for immunoblot analysis of matrix metalloprotcinase-3 
(MMP-3) and A disintegrin and metalloproteinase with thrombospondin motifs 5 (ADAMTS-5). Cell 
layers were extracted for analysis of nuclear factor kB (NF-kB) and inducible nitric oxide synthases 
(iNOS). The results revealed that lysates, Fn-f, and IL-1/3 induced NF-kB phosphorylation (Figure [T]A) . 
Fn-f, and IL-1/3 but not lysates also induced iNOS, MMP-3 and ADAMTS-5 expression (Figure[l^,C,D). 
In the case of NF- B, EPO treatment blocked activation by lysates, but not by Fn-f or IL-1 . EPO also 
dramatically reduced iNOS, MMP-3, and ADAMTS-5 expression in response to Fn-f and IL-1/3. These 
results underscore the ability of EPO to counter injury effects on multiple pathways. 

Figure [2] illustrates the balance between pro-inflammatory (e.g. TNF-a) and anti-inflammatory (e.g. 
EPO) cytokines during the typical injury response. Damage due to chemical or mechanical stress initiates 
production of alarmins, such as damage associated molecular pattern molecules (DAMPs) [5]-[7]. These 
alarmins trigger production of pro-inflammatory cytokines such as TNF-a by cells near the initial injury 
site (see (b) in figure [2| Thus there is a penumbra ((c) in figure [2| of "sick" cells formed at the 
boundary of the initial injury as a result of inflammation. While still viable, these cells are at risk of 
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Figure 1: Effects of EPO on chondrocyte responses to alarmins and IL-1/3. Immunoblot 
analyses are shown for phospho- and total nuclear factor kappa-B (NF-kB), a pro-inflammatory tran- 
scription factor (A), inducible nitric oxide synthase (iNOS), a catabolic signaling factor (B), matrix 
metalloproteinase-3 (MMP-3) (C), and a disintegrin and metalloproteinase domain with thrombospondin 
repeats (ADAMTS-5) (D), which cause cartilage matrix degeneration. Primary cultures of bovine chon- 
drocytes (1 X 106 cells/well) were treated with 10 lU/ml EPO for 24 hours, then challenged for an 
additional 24 hours with the indicated concentrations of purified 29 kDa fibronectin fragment, a known 
extracellular matrix DAMP, or interleukin-lb (IL-lb), a pro-inflammatory cytokine, or cell lysates con- 
taining multiple cellular DAMPS. EPO treatment suppressed catabolic responses induced by all three 
treatments. 



dying themselves due to the presence of the pro-inflammatory cytokines, resulting in the spread of injury 
as shown in (d)-(e) in figure [2j However, cells far enough from the injury are able to produce anti- 
inflammatory cytokines such as EPO to check the spread of injury and promote healing. This is shown 
(d*)-(e*) in figure[2] It is also suggested in [5' that TNF-a acts in such a way to limit EPO signaling. Thus 
there is a critical balance between pro-infiammatory and anti-inflammatory cytokines which determines 
the spreading properties of an articular cartilage lesion associated with injury. Moreover, the authors 
in [5] suggest that introduction of exogeneous EPO at an early time after an initial injury could be an 
effective therapy for minimizing the amount of collateral damage to cartilage that occurs as a result of 
inflammation. 

The goal of this work is to study via mathematical modeling the EPO/TNF-a interaction described 
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Figure 2: Injury response. When a wound occurs DAMPs promote production of TNF-a (b), an 
initially viable penumbra (yellow in (c)) of cells becomes EPOR active. These cells are in danger of death 
which can result in the spread of injury as seen on the left in (d) and (e). However, EPO can counteract 
the effects of TNF-a and prevent spread of the lesion (d*) and (e*). 

above. Herein we develop a model capable of simulating the scenario illustrated in figure [2j Moreover, 
the model is developed as a potential for studying the use EPO or other anti-inflammatory cytokines as 
a therapy for lesion abatement. The type of model, and in particular the one described below, in this 
article is of value as it provides a link between observed phenomena and the mechanisms involved in 
driving that phenomena. Models such as those presented below can also help to generate hypotheses and 
suggest experiments or fundamental quantities to be measured. More practically, mathematical modeling 
such as we use in this paper could provide quick and inexpensive initial screening of potential therapies. 



Models 

In this section we describe the mathematical model developed to represent the biological interactions of 
chondrocytes and cytokines described in [s]. We aim for a minimal model based on known mechanisms 
considered to be the dominant factors in articular cartilage lesion abatement. By a minimal model 
we mean one in which the removal of any component results in behavior that is inconsistent with the 
typical injury responses in cartilage as discussed in [5]. We note that the model described below can 
be expanded or modified to include further interactions as new results from experiment and observation 
dictate. In particular, further chemical pathways such as interleukin-1 j3 IL1-/3 can easily be incorporated 
if necessary. The chemical species included below are chosen for their functionality, that is their action on 
chondrocytes during injury response. Hence they can be replaced with any other chemical species whose 
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Figure 3: States of chondrocytes during stereotypical injury response. Healthy cells (C) become 
sick. Typically a healthy cell will first become catabolic {St) and either enter apoptosis (Da), or become 
EPOR active (Sa)- EPOR active cells may be saved by EPO to become healthy once again, or die 
(apoptosis). 



effects on chondrocytes are functionally the same. The chondrocyte cell states described below represent 
the biological actions of the cells in response to signaling by a particular cytokine during the typical 
injury response. These actions are assumed to be analogous to those of other cell types, for which aspects 
of innate immune response, such as local inflammation considered in this paper, are well established. 

We consider a population of chondrocytes fixed in matrix, which is typical of articular cartilage. Sub- 
populations of chondrocytes are assumed to exist in different states corresponding to the chemical signals 
being received by the cells during injury response. Figure |3] shows the states in which subpopulations of 
chondrocytes may exist. We refer to the normal state of a subpopulation of chondrocytes as the healthy 
state. We denote by C the population density (cells per unit area) of healthy chondrocytes at a given time 
and location. As a result of inflammation and injury healthy chondrocytes can enter into a "sick" state. 
Cells in this state are at risk death (via apoptosis) unless their signaling by TNF-a is somehow limited. 
We consider two subpopulations of cells in the sick state. We denote by St the population density of cells 
in the "catabolic" state. Catabolic cells are chondrocytes that have been signaled by alarmins and are 
capable of synthesizing TNF-a and other cytokines associated with inflammation. Healthy cells signaled 
by DAMPs or TNF-a enter into the catabolic state and begin to synthesize TNF-a and produce reactive 
oxygen species (ROS). Catabolic cells that are signaled by TNF-a express a receptor (EPOR) for EPO 
and make up the subpopulation of sick cells we refer to as EPOR active. It should be noted that there 
is a time delay of 8-12 hours before a cell expresses the EPO receptor after being signaled to become 
EPOR active [Sj. We denote the population density of EPOR active cells by Sa- Since EPOR active cells 
express a receptor for EPO, they may switch back to the healthy state if signaled by EPO. However, as 
discussed in [s] TNF-a limits production of EPO. Thus there is a balance between EPO and TNF-a that 
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Figure 4: Signaling involved in cartilage injury response. The occurrence of an injury begins 
a sequence of chemical productions that promote the inflammatory response. Lysing (necrotic) cells 
give off DAMPs (M) resulting in a population of catabolic cells that produce TNF-a (F) and ROS 
{R). ROS promotes production of EPO (P) by healthy cells while TNF-a acts to block EPO signaling. 
Furthermore, TNF-a degrades extracellular matrix providing another source of DAMPs. The cells switch 
state according to the signaling by chemicals produced in response to injury and inflammation. DAMPs 
and TNF-a drive the C ^ St transition, TNF-a drives the St — >■ Sa transition, while TNF-a and 
DAMPs drive the St Da transition. Finally, EPO drives the Sa ^ C transition. This figure 
illustrates the assumptions of the mathematical model ((|3| below). Not shown: DAMPs given off as a 
result of matrix degradation also influences the switch from healthy to catabolic. 

determines the spreading behavior of cartilage lesions. The catabolic and EPOR active cells together 
make up the population of cells forming the penumbra as illustrated in figure [2] We also consider a 
"dead" state for subpopulations of chondrocytes. This includes necrotic cells ZJjv and apoptotic cells Da- 
We note that for the purposes considered herein that apoptotic cells do not feed back into the system 
and thus are not explicitly represented in the mathematical model. We assume chondrocytes become 
necrotic only at the time of an initial injury. Necrotic cells release alarmins such as DAMPs that initiate 
the injury response. Either catabolic or EPOR active cells become apoptotic if signaled by inflammatory 
cytokines. 

Figure |4] details the chemical signaling and the switching of chondrocyte cell states represented in the 
mathematical model presented in the following. An initial injury creates a population of necrotic cells 
which release alarmins (such as DAMPs) [6|[7j . We denote by M the concentration of DAMPs at a given 
time and location. The DAMPs signal healthy cells near the injury to enter the catabolic state resulting 
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in the production of TNF-a. We denote the concentration of TNF-a at a given time and location by F. 
The inflammatory cytokine TNF-a has several effects on the system: It 

1. feeds back to continue to switch healthy cells into the catabolic state, 

2. causes catabolic cells to enter the EPOR active state [s], 

3. influences apoptosis of catabolic and EPOR active cells, 

4. degrades extracellular matrix (denoted by U) which results in increased concentrations of DAMPs, 

5. has a limiting effect on production of EPO [5]. 

Catabolic cells also produce reactive oxygen species (ROS) which influences the production of EPO by 
healthy cells. We denote the concentration of ROS at a given time and location by R. There is a time 
delay of 20-24 hours before a healthy cell signaled by ROS will begin to produce EPO. 

In developing a mathematical model to represent the scenario described above we assume that the 
chemicals diffuse throughout a domain. This diffusion is essential in determining the spatial behavior, 
i.e. the expansion or abatement of the lesion. Since the chondrocytes are fixed in the matrix we do 
not consider cell motility. First the mathematical model contains four equations (one for each chemical 
species) describing the dynamics of the chemical concentrations. These equations are each of the form: 

change in concentration of particular chemical 
= diffusion of that chemical 

+ production of that chemical by cells in appropriate states 

— natural decay of that chemical. (1) 

Next the model consists of four equations for the population densities of cells in the healthy, catabolic, 
EPOR active, and necrotic states. Each of these equations are of the form: 

change in population density of particular cell state 

= population density of cells switching into that state due to appropriate signaling 

— population density of cells switching out of that state due to appropriate signaling. (2) 
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Finally there is an equation corresponding to the degradation of extracellular matrix by TNF-a which 
feeds back into the production of DAMPs. 

The model equations for the concentrations of the chemical species are: 

dtR^V ■ {DrVR)-5rR + (jrSt, (3a) 

dtM =V • {DmVM) - 5mM + <JmDn + SuU—^, (3b) 

dtF ■ {DfVF) -SpF + <7fSt. (3c) 

a^P =V . (DrVP) - SfP + aFdt - r.)j£l^^^^. (3d) 

The equation describing matrix degradation by TNF-a is: 

d,U = -6uU—^. (3e) 



We note that the right hand side of this equation appears in ( 3b I as part of the production of DAMPs 
as we have assumed that degraded matrix releases alarmins. 
The equations for the switching of cell states are: 

P K'f P 

-as,— - hC^-^^H^P - P.) - hC—HiP - P.), (3f) 

F M 
-"^"-Xf + FXm + M' 

dtDN ^ - IIDmDn- (3i) 



In (kot and (3gl the function H{-) is given by 



1 if s < 0, 

H{s) = { (4) 
if s > 0. 

The constant P^ represents a critical level of EPO above which the effects of DAMPs and TNF-a on 
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healthy cells is limited. 

Based on the diagram in figure |4] and the assumptions underlying that diagram, one could assume 
that the function H{-) is identically one, i.e. H = 1. However, if we take H = 1 then figure |4] implies, 
and computational results confirm, that there is infinite feedback into the system by alarmins, TNF- 
a, and catabolic cells. As noted in [s] "the pro-inflammatory arm of the injury response is inherently 
self-amplifying" . This does not allow for lesion abatement or limitation of secondary pro-inflammatory 
cytokine induced injury described in [Hj. This suggests that signaling by anti-inflammatory cytokines such 
as EPO not only promote the switch from the EPOR active state to the healthy state but also influences 
the response of healthy cells to alarmins and pro-inflammatory cytokines. Thus the mathematical model 
suggests that in some way the anti-inflammatory cytokines limit the switch from the healthy state to the 
catabolic state so that there is not an infinite feedback into the system by alarmins, TNF-a, and catabolic 
cells. This is consistent with the observation of Brines and Cerami on the antagonistic relationship 
between TNF-a and EPO, that "each is capable of suppressing the biological activity of the other" 
In the model equations ([3| we take the function H{-) as in for convenience. However, there may be 
a more appropriate form for the function H{-) that must be determined through experiment to discover 
the nature of the effects that sufficient concentrations of EPO or other anti-inflammatory cytokines have 
on healthy chondrocytes. 

We note some features of the mathematical model ID 

1. we have incorporated the time delays for activation of the EPO receptor and synthesis of EPO 



(3dM3g),(3h) 



2. it requires a concentration of TNF-o; and DAMPs together for apoptosis of catabolic cells (see (3g)) 



3. after some time necrotic cells decay to an inert state (see ([3i])) so that there is not a continuous 
production of DAMPs for all time from the initial injury. This also corresponds to the loss of 
cartilage such as is sometimes associated with osteoarthritis. 

The speciflc functional forms appearing in the model system ([3]) have been chosen to capture the critical 
thresholds and represent the function shapes qualitatively. 
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Results and Discussion 

Using the mathematical model ([3| we can simulate the two scenarios depicted in figure [2] That IS, 
we can choose parameter values that correspond first to the case where infiammation is not hmited by 
the presence of anti-infiammatory cytokines (d) and (e) in figure [2] and second to the case where anti- 
infiammatory cytokines hmit the development of secondary inflammation induced injury (d*) and (e*) 
in figure [2j For the computations and simulations we assume a two dimensional dimensional domain. 
Moreover, as depicted in figure [2] we assume circular symmetry so that changes only occur in the radial 
direction. In order to simulate the injury response we obtain numerical approximations to the system 
of partial differential equations ([3|. This is done by discretizing the system in the spatial terms using a 
standard finite difference scheme appropriate for circular symmetry. This results in a system of delay- 
differential equations which are solved in MATLAB using the dde23 code for solving delay-differential 
equations. For details on the methods and software for solving delay-differential equations see [8-10 . 



First, using the mathematical model ([s]) we simulate the evolution of chondrocyte population density 
(per area, and assuming circular symmetry, as a function of radius) over a period of ten days after an 
initial cartilage injury. This corresponds to the cascade of events illustrated in the left hand arrows of 
figure[2] The initial conditions correspond to the presence of an initial injury, i.e. a population of necrotic 
cells occupying a circle of radius 0.25 mm. 

Following the initial injury there is a quick and steady increase in the population density of catabolic 
and EPOR active cells and decline in the population density of healthy cells (figures [5(b)p(i)| ). After ten 



days there has formed a thick penumbra of catabolic cells which are beginning to die off (via apoptosis. 



see figure 5(i) ) showing that without the presence of EPO there is unlimited secondary pro-inflammatory 
cytokine-induced injury and unabated spread of the cartilage lesion. This corresponds qualitatively well 
with the expected results illustrated in figure [2] (d),(e). 

Figure [g] shows the time evolution of chondrocyte population density (per area, and assuming circular 
symmetry, as a function of radius) over a period of ten days after an initial cartilage injury with the 
antagonistic effects of EPO on TNF-a. The initial conditions are the same as previously described for 
the simulations without EPO. The dynamics of the chondrocyte cell states population densities are the 
same over the first day or so as in the case with no EPO. This is due to the time delays for cells to 
become EPO producing and EPOR active. However, we see by day 3 or so (see figure [6(d)] and compare 
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with 5(d) ) there is a significant difference in the population dynamics compared with that in the no EPO 
case. We begin to see a reversion of cells in the EPOR active state back to the healthy state. While 
there is still some amount of secondary injury due to the inflammatory process and presence of pro- 
inflammatory cytokines, there is a fixed radius (see figures |6(f)|6(i)" ) beyond which the penumbra, and 
hence the lesion and secondary pro-infiammatory cytokine induced injury, cannot expand. Thus there is 
lesion abatement by the action of the anit-inflammatory cytokine EPO. This corresponds qualitatively 
well with the expected results illustrated in figurc[2] (d*),(e*). 

Based on the biological assumptions summarized in the diagram shown in figure[2]we have developed a 
mathematical representation of some of the principal features of the injury response in articular cartilage. 
In particular we have captured, qualitatively, secondary pro-inflammatory cytokine induced injury and 
lesion based on the balance between TNF-a and EPO as described in Brines and Cermani [5]. This mod- 
eling effort suggests that anti-inflammatory cytokines may act in such a way as to influence the response 
of chondrocytes to alarmins and pro-inflammatory cytokines. In particular sufficient concentrations of 
anti-inflammatory cytokines may limit healthy chondrocytes conversion to the catabolic state. Provided 
dosages, dosage responses, delivery method, etc. the model ([s]) can be applied to study the effectiveness 
of the introduction of exogenous EPO as a potential therapy for limiting secondary pro-inflammatory 
cytokine induced injury and promoting healing in articular cartilage. We note however, that due to the 
time delays for cells to become EPOR active, there is little or no response to EPO at early times after an 
injury has occurred. Thus there is a window of time during which introduction of tissue-protective EPO 
derivatives such as discussed in |5 are ineffective as therapies. This is evidenced by comparing flgures 
5(a)||5(b) with figures 6(a)|6(b) in which there is little or no difference between the case with EPO and 



the case without. The model presented herein also has potential applications for studying other issues 
related to post-traumatic stress in articular cartilage. By coupling this model with mechanical models, or 
including geometric features of cartilage one may be able develop a more complete theoretical description 
of aspects of cartilage injury of relevance in the biomedical sciences. 
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Figure 5: Healthy and Penumbra cells without EPO. This figure shows the evolution of the healthy 
and penumbra (comprised of both catabolic and EPOR active) cell population densities over a period 
of 10 days. The results shown here correspond to a situation in which there is no production of or 
response to EPO. We see here that the healthy cell population is decimated. In accord with figure |4] and 
the corresponding mathematical model ([s]) this is due to the uninhibited influence of pro-inflammatory 
cytokines. 
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Figure 6: Healthy and Penumbra cells with EPO. This figure shows the evolution of the healthy 
and penumbra cell population densities over a period of 10 days. Initially there is the development of 



a population of penumbra cells (6(c)) due to the time lags for EPO production and EPOR activation. 



However, once healthy cells at the edge of the penumbra begin to produce EPO and penumbra cells 
become EPOR active there is an increase in the population density of healthy cells near the injury and 



there ceases to be expansion of the penumbra (figures 6(d) - 6(i) ) 



14 



References 

1. Buckwalter J, Mow V, Ratcliffe A (1994) Restoration of injured of degenerated articular cartilage. 
J Am Acad Orthop Surg 2: 192-201. 

2. Freedman K, Coleman S, Olenac C, Cole B (2002) The biology of articular cartilage injury and 
the microfracture technique for the treatment of articular cartilage lesions. Seni Arthroplasty 13: 
202-209. 

3. Goodwin W, McCabe D, Sauter E, Reese E, Walter M, et al. (2010) Rotenone prevents impact- 
induced chondrocyte death. J Orthop Res 28: 1057-1063. 

4. Ramakrishnan P, Hecht B, Pedersen D, Lavery M, Maynard J, et al. (2010) Oxidant conditioning 
protects cartilage from mechanically induced damage. J Orthop Res 28: 914-920. 

5. Brines M, Cerami A (2008) Erythropoietin-mediated tissue protection: reducing collateral damage 
from the primary injury response. J Intern Med 264: 405-432. 

6. Harris HE, Raucci A (2006) Alarming news about danger. EMBO Reports 7: 774-778. 

7. Bianchi ME (2006) Damps, pamps and alarmins: all we need to know about danger. J Leukoc Biol 
81: 1-5. 

8. Shampine LF, Reichelt MW (1997) The matlab ode suite. SIAM J Sci Comput 18: 1-22. 

9. Shampine LF, Thompson S (2001) Solving DDEs in MATLAB. Appl Numer Math 37: 441-458. 
10. Shampine LF, Thompson S (2009) Numerical solution of delay differential equations : 245-271. 



